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Abstract 

Surface-induced profiles of both nematic and smectic order parameters in a nematic liquid crystal, 
ranging from an orienting substrate to "infinity" , were evaluated numerically on base of an extended 
Landau theory. In order to obtain a smooth behavior of the solutions at "infinity" a boundary 
energy functional was derived by linearizing the Landau energy around its equilibrium solutions. 
We find that the intrinsic wave number of the smectic structure, which plays the role of a coupling 
between nematic and smectic order, strongly influences the director reorientation. Whereas the 
smectic order is rapidly decaying when moving away from the surface, the uniaxial nematic order 
parameter shows an oscillatory behavior close to the substrate, accompanied by a non-zero local 
biaxiality. 
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I. INTRODUCTION 



The structure of uniaxial nematic liquid crystals is strongly influenced by the presence 
of orienting surfaces [1]. In addition to the usual elastic distortions in the nematic bulk, 
which are well described by the Oseen-Zocher-Frank energy, strong deformations can occur 
close to the surface substrate [2, 3]. These effects are often accompanied by a non-zero local 
biaxiality of the orientational order [4]. The nematic orientation close to a confining sub- 
strate could be detected experimentally using second-harmonic generation techniques [5]. In 
addition, by X-ray studies the surface has also been proven to induce a layered structure, 
i.e., smectic order appears close to the surface which decays when moving away from the 
substrate into the nematic bulk [6]. In spite of considerable effort, both in theory [7, 8, 9, 10] 
and computer simulation [11, 12, 13, 14], the full complexity of surface-induced structural 
changes in nematics is far from being understood. Whereas Monte Carlo or molecular dy- 
namics simulations approach the problem on the molecular level, in our contribution we 
take a phenomenological viewpoint. To this aim we consider an extended Landau theory, 
comparable to Skacej et al. [8]. In addition to that paper, we investigate not only uniaxial 
order, instead, both the full alignment tensor for nematic order and, mainly, the amplitude 
of the smectic layering is taken into account. In particular, the influence of the coupling 
between smectic and nematic order on the order parameter profiles, obtained from numerical 
relaxation, is investigated in detail. Whereas the boundary conditions at the surface sub- 
strate are fixed, an additional boundary energy is derived, in order to guarantee a smooth 
behavior of the profiles at "infinity" (far from the surface) where the volume equilibrium 
values of the order parameters should be reached. 

The organization of the article is as follows. In Section II the Landau theory used in our 
calculations is introduced. (We focus on the derivation of the additional boundary energy at 
"infinity" in the Appendix.) Section III indicates the numerical relaxation method employed 
and presents selected results for order parameter profiles. Finally, Section IV contains some 
concluding remarks. 
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II. EXTENDED LANDAU THEORY FOR SURFACE-INDUCED EFFECTS 



The geometry of our system is the semi-infinite space (z > 0), confined by a substrate 
surface at z = and "infinity" (z = oo) which, in numerical practice, means a large distance 
from the surface. Due to the infinite extension of the system in x and y direction and 
the absence of any lateral structure of the surface, we can reduce the problem to a one- 
dimensional geometry, i.e., all quantities only depend on the distance z from the surface. In 
order to be able to investigate both positional and orientational order, we need two different 
order parameters. The smectic order parameter usually is a complex number, p = ipe lx , 
whose phase x accounts for local layer deformations. We assume perfect layering at the 
surface (p(z = 0) = 1) and, therefore, we are left only with the amplitude ip(z) of the 
layering which is a real quantity indicating the degree of smectic order. The nematic order 
parameter is a second-rank traceless and symmetric tensor. Without loss of generality we 
choose its parametrization as 



(i) 



Qxy\%) Qyy\%) Qyz(^) 
^ Qxzi^z) Qyz\Z) Qxx\Z) Qyy\Z) j 

Therefore, there are six scalar functions whose profiles (z dependence) have to be deter- 
mined. These profiles are found from a numerical minimization of an energy, supplied with 
appropriate boundary conditions at z = and z = oo. 



A. Bulk energy functional 

The bulk energy functional is chosen according to an extended Landau theory. It consists 
of a smectic and a nematic contribution. The smectic energy contains a volume and an elastic 
contribution [15], 

^smec = -r|p| 2 + \ P \ 4 + -k\(X7 -iq 5n ± )p\ 2 . (2) 

2 2 

Due to the continuity of the smectic-nematic phase transition, the volume smectic energy is 
an expansion into even powers of the smectic order parameter p. This phase transition occurs 
at t = 0, where r is a reduced temperature. The last expression in (2) is the elastic energy 
due to the gradients of the layer amplitude (with smectic elastic constant k). It incorporates 
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a coupling to the nematic order, based on local U(l) gauge invariance, due to the nature of 
the smectic order parameter as a complex number [15]. Namely, local changes of the smectic 
order are accompanied by transversal fluctuations of the director field n, perpendicular to 
the layer normal. In our simplified geometry these fluctuations are always in the x-y plane. 
The coupling strength is given by the intrinsic wave number q of the smectic layering. After 
inserting the nematic tensor order parameter (1) and reducing the smectic order parameter 
to its amplitude, the smectic energy becomes a functional dependent on ip(z), ip' (prime 
denoting derivative with respect to z) and Qij(z), 



1 i.2 , „,4 1 



xx + Qyy "I" S vo i) ^] • (3) 



Here, S vo \ is the scalar order parameter that minimizes the volume nematic energy (see next 
subsection) . 

The nematic energy functional also consists of a volume and an elastic part, 

3~ nem — t Qij Qji V^6 Qij Qjk Qki H"~ (Qij Qji) (4) 

4 

+- (diQjk) (diQjk) + - k 2 i (diQik) (djQ jk ) + - fc 31 (d k Qij) (djQ ik ) 

4 4 4 

~ t {Qxx Qxy Qxz Qyy Qyz ~\~ Qxx Qyy) (5) 



3 \/6 (Qxx Qxy ~t~ Qyy Qxy Qxx Qyy Qyy Qxx 
Qxx Qyz Qyy Qxz 2 Qxy Qxz Qyz^) 

(Qxx Qxy Qxz Qyy Qyz Qxx Qyy) 

+ ~ (Qxx + Qxy + Qxz + Qyy + Qyz + Qxx Qyy) 



+- (k 21 + k 3l ) ■ (Q' xx 2 + Q' yy 2 + Q' xz 2 + Q' yz 2 + 2Q' xx Q'J. 



Unlike the smectic energy, the volume part of (4) contains a third-order expression, to 
describe the discontinuous isotropic-nematic phase transition, which, in our parametrization. 
occurs at t = §. (f is again a reduced temperature.) For the elastic part of (4) there are 
three independent deformation modes, similar to the Oseen-Zocher-Frank theory for elastic 
distortions of the director. k 2 \ and k 3 i denote the ratios of elastic constants for the respective 
deformation modes. 
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B. Boundary conditions 



At the substrate surface the values for the order parameters are fixed. We assume ideal 
smectic and uniaxial nematic order. That means the smectic amplitude is one, and the 
alignment tensor is completely determined by the uniaxial scalar order parameter 5 = 1 
and the fixed director surface tilt angle 6 sur f (measured in the x-z plane, from the z axis). 
This results in the Dirichlet boundary conditions, 



#3 = 0) = 1 , 

Q(z = 0) = S (n®n--l) 
/ 



sin 6 surf - 3 


sin GL, r f cos 



surf 



sin 6 surf cos surf 

-I o 

cos 2 surf - \ 



\ 



J 



(6) 
(7) 

(8) 



At infinity, the boundary conditions are not of Dirichlet-type. Instead, we have to guar- 
antee smooth profiles (zero slope) for the order parameters which should reach those values 
that minimize the volume parts of the smectic and nematic energy (3) and (5). To this aim 
we insert the uniaxial form of the alignment tensor (7) into (5). A direct minimization yields 
the temperature dependence of the volume order parameters, 



VVoi(t) = 
S vo i(t) = 



2 V ~ T 



3 V6 + V54 - 48 1 
16 



(9) 
(10) 



The tilt angle to be reached at infinity, 6oo, is unknown a priori. (Previous molecular 
dynamics simulations indicate a tilt angle of ©oo = at infinity [12].) We, therefore, take 
6oo as a free parameter in our calculations. Its actual value is determined by performing a 
series of simulations for any fixed set of the remaining simulation parameters. Monitoring 
©oo versus the energy E cqui of the equilibrated configurations for this series yields the tilt 
angle at infinity, which corresponds to the minimum of the function Q 00 (E equ i). Following 
this procedure, we are able to determine Goo in dependence of the reduced temperature t 
and the smectic wave number q Q . 

Next we have to establish a boundary energy functional whose minimization leads to 
the desired surface profiles at infinity. Based upon the procedure introduced by Galatola et 



al. [16] the central idea for finding the boundary energy functional is a linearization of the 
bulk equations around the volume order parameters calculated above. We provide all details 
of this derivation, which is an important ingredient of our method, in the Appendix. Here 
we merely quote the result, 

^00 = \U (Q xx 2 + Qyy 2 ) + - 3 (U + V) S vol (Q xx + Q yy ) (11) 
SVol sin ©oo ' (U Qxx ~\~ V Qyy) V Qxx Qyy H • Qxy 

+ 1 - (2 + k 2 i + hi) t ■ (Qxz 2 - 2 SVoi sin 0^ cos ©^ • Q xz + Q yz 2 ) 

+ - y/KT ■ 1p 2 , 
2 

where U and V are terms which depend on the reduced temperatures r and t, the ne- 
matic elastic constants k 2 i, A: 31 and the volume order parameter S vo \(t). (For the explicit 
expressions see Appendix A.) 

Now all order parameter profiles can be obtained from a minimization of the total energy 

POO 

F= (^ + ^ 00 5(oo))dz , (12) 
Jo 

with the Dirichlet boundary conditions (6) and (8) valid at z — 0. 



III. ORDER PARAMETER PROFILES 

The minimization of the total energy (12) was performed numerically, employing a stan- 
dard Newton- Gaufi-Seidel technique which, in our case, formally corresponds to a one- 
dimensional version of the Finite Element method. First the "infinite" distance from the 
surface was replaced by a large, finite value z max = 100. The range < z < z max was 
discretized in N = 1000 intervals. The bulk and boundary energy functionals (3), (5), (11) 
were evaluated on these intervals, the derivatives with respect to z being replaced by finite 
differences. The values of the order parameters at z = were fixed according to (6) and 
(8). For the initial configuration we assumed linear profiles for all quantities ip and Qij 
on the interval < z < z max , by interpolating between their surface and volume values. 
An iterative procedure was then performed on each grid point. All order parameters were 
corrected according to the Newton-GauB-Seidel prescription, 

xncw = x ° ld - ' (x = ^ q ^ ' (i3) 
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where the functional derivatives in (13) were evaluated by numerical differentiation. The 
relaxation was terminated when the relative change of the total energy was less than 10~ 6 
which corresponded to some thousand relaxation steps. 

For further discussion the nematic tensor order parameter will be analyzed in terms of 
its eigenvalues and eigenvectors. More specifically, we will plot the tilt angle of the main 
director, measured from the z axis and two scalar order parameters. The latter ones measure 
the degree of uniaxial and biaxial order, respectively. (We checked that the director twist 
angle stays constant, due to the surface anchoring in the x-z plane.) 

The equilibrium profiles were evaluated for different values of the reduced temperature 
t = 0...1 and the intrinsic smectic wave number g = 0.3. ..0.8. All other parameters were 
fixed, except for the tilt angle at infinity Ooo, which was found for each set (t, go) from 
the additional minimization procedure indicated in Section II. B. Thereby, could be 
determined up to a maximum error of ±1°. Regarding the remaining parameters, the elastic 
constants were chosen as fc 2 i — k 3 ± — 1 and k — 5 which accounts for the fact that layer 
distortions should contain a higher elastic energy than deformations of orientational order. 
The reduced temperature was r = 0.1 which corresponds to a nematic state point, slightly 
above the smectic-nematic phase transition. Finally, the director at the surface was anchored 
at a tilt angle of 6 sur f = 60°. 

Let us first discuss the behavior of the tilt angle at infinity 6oo in dependence on the 
reduced temperature t and the smectic wave number q . The respective results are displayed 
in Table I. Fixing the surface tilt angle at © sur f = 60°, in almost all cases we find a strong 
reorientation towards the surface normal which increases with the smectic wave number. 
The influence of the smectic wave number on the director reorientation can be understood 
from the particular form of the coupling energy (2). In order to minimize this coupling the 
transversal director components should be small in those regions where the smectic order 
parameter is significantly non-zero, i.e., close to the surface. Therefore, it is obvious that for 
increasing wave number the reorientation of the tilt angle towards the direction of the surface 
normal becomes more pronounced. Unlike for the smectic wave number, the dependence of 
Ooo on the reduced temperature t is fairly small. Only for low values of both t and go the 
tilt angle deformation is reduced. Apparently, at t — 0.0 the orientational fluctuations are 
still too small to enhance the director reorientation for small wave numbers. 

Figures 1-3 correspond to a reduced temperature of t — 0.0 which means a nematic state 
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point far away from the nematic-isotropic phase transition. The smectic order parameter 
profile is given in Fig. 1 for wave numbers go = 0.3, 0.6, 0.8, corresponding to a layer spacing 
of 20.9, 12.6, 7.8. Obviously, the smectic structure is rapidly decaying from its maximum 
value of one when moving away from the surface. At a distance of 20 the curves essentially 
have reached their asymptotic value of zero. This loss of smectic order is almost independent 
of the wave number. 

As shown in Figure 2 the uniaxial nematic order parameter is also approaching its asymp- 
totic value within a distance of 20. The value reached is 0.92, which precisely corresponds 
to the volume value of the scalar order parameter at temperature t = 0.0, according to (10). 
Remarkably, the uniaxial order parameter is not decaying monotonically, instead there is an 
oscillatory behavior. For large wave numbers (q > 0.6) it even decreases below its volume 
value. The non-monotonic behavior of uniaxial order is accompanied by the occurence of 
a non-zero biaxial order parameter. Both suppression of uniaxial order and increased bi- 
axiality close to the surface have also been observed in computer simulations based on the 
molecular Gay-Berne model [12]. 

Unlike the scalar order parameters, the profile of the director tilt angle strongly depends 
on the intrinsic smectic wave number (Figure 3). Whereas for q = 0.8 there is again a strong 
change within a comparably short distance, for low wave numbers (q < 0.6) the director 
reorientation is much weaker. The behavior of the tilt angle is also changing qualitatively 
with the wave number. E.g., for q = 0.8 the tilt angle profile becomes non-monotonous, 
taking intermediate values that are closer to the homeotropic orientation than the tilt angle 
6oo finally reached. As revealed from Figure 4, an even more drastic change occurs for the 
reduced temperature t = 1.0, which is just below the nematic-isotropic phase transition. For 
low wave numbers the tilt angle profiles show a local maximum at a distance of around 20, 
before decaying towards the volume value. 

IV. REMARKS 

1. Summarizing our work, we have numerically analyzed the surface-induced profiles of 
smectic and nematic order as well as director orientation. Whereas the order param- 
eters are always strongly changing in a thin layer close to the surface, the tilt angle 
reorientation is mainly dependent on the intrinsic smectic wave number of the liquid 
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crystal which, in our model, acts as a coupling parameter between nematic and smectic 
order. It seems, however, that the case of high wave number (q = 0.8) is the most re- 
alistic, considering the experimental observation of a strong reorientation close to the 
surface [5]. In addition, surface-induced biaxiality and suppression of uniaxial order, 
previously detected in experiment and molecular simulations, could also be confirmed 
within the frame of our model. 

2. The asymptotic behavior of the profiles at infinity was reproduced in numerics by 
deriving a boundary energy functional. To this aim, referring to the method of Ref. [16] 
we linearized the bulk Euler-Lagrange equations around their volume solutions. When 
the boundary energy is included into the numerical minimization procedure, the order 
parameter profiles at infinity smoothly reach their volume values. Although in our 
case this technique was used in a one-dimensional geometry, it could be extended in a 
straightforward way to more complex situations. For instance, a laterally structured 
surface, which is fairly common in novel developments of display technique, will already 
break the symmetry of our example. However, the treatment of open boundaries by an 
additional energy functional derived from a linearization of the bulk energy is a very 
general concept which could be used for any numerical study dealing with "infinitely" 
extended systems. 

APPENDIX A: DERIVATION OF THE BOUNDARY ENERGY FUNCTIONAL 

Before starting to derive the boundary energy functional at infinity, we first remind of 
the desired values for the order parameters. At infinity the smectic order is completely lost 
(ip(z = oo) = 0). The nematic order parameter is again uniaxial, but with the volume scalar 
order parameter S vo \(t), dependent on temperature, and tilt angle Goo, 



i]){z = oo) = 



(Al) 




Svo\(t) (sin 2 G 







S vo \(t) sin Gqo sinG 



\ 



oo 



CO 








?vol(0 

S vol (t) (cos 2 G, 
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oo 



oo 
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We start with the bulk equations which are obtained from the bulk energy functional 
J~ = ^smec + ^nem, (3) and (5), by variational calculus as the corresponding Euler-Lagrange 
equations, 



dT d dT 



(A3) 



ax d^ = °- = 

In order to obtain a smooth behavior of the profiles at infinity, the Euler-Lagrange equations 
(A3) are linearized around the volume order parameters. To that aim we calculate the partial 
derivatives of the right-hand sides at infinity which, in turn, are obtained from the quadratic 
part of the free energy (3) and (5), 
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The linearized Euler-Lagrange equations now become explicitly 

(3 + fc 21 + k 31 )t 



Qxi 



3 + 2(hi + hi) K-^oi (^©oo-^)) 
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(Qxx ~ S vo i (sin 2 ©oo - ^)) 
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(Qxz - S vo i sin 6oo sin 0^) , 
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K 



t , 



1 

2* 



t. 



(A4) 
(A5) 
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Among the equations (A7)-(A12), the first two are coupled. In order to find their solu- 
tions we note that they can be written in matrix form, 




A B\ I Q xx - SVoi (sin 2 Ooo - §) 

B A J \ Qyy "I" \ »Svol 



(A13) 



' yy 

Now we change the variables to the deviations of Q xx and Q yy from the volume solutions, 

(A14) 

By determining the eigenvalues and eigenvectors of the coefficient matrix, it can be expressed 
by its diagonalized form and the corresponding orthogonal transformations, 






(A15) 



SQ'L 

Wyy 

Now the system can be decoupled by a similarity transformation. Keeping only the decaying 
modes, we arrive at the solution 

/ 1 1 \ / e ~VA=Bz o 

■11 e -^+B-' 




1 

V2 
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(A16) 



The asymptotic Cauchy boundary conditions are obtained by differentiating the solutions 
by z, taken at infinity, 



5Q' XX ] = _ J_ [ 1 1 ] . f \ j_ / 1 -1 

The explicit boundary conditions for Q xx and Q yy then are 




(A17) 



Q' xx (z = 00) = - 1 - (VX+B + VA^j ■ (Q xx - S vol (sin 2 ©00 - ^)) (A18) 
- 1 - (VA + B - VA^j ■ (Q yy + l - S vol ) , 

Q' yy {z = 00) = -- (VA + B - VA^B^j ■ (Q xx - S vol (sin 2 ©00 - ^)) (A19) 
- 1 - (VATB + VA^B^j ■ [Qyy + l - sjj . 
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The remaining linearized Euler-Lagrange equations are already decoupled which immediately 
yields the corresponding Cauchy boundary conditions at infinity, 

Q' xy (z = oo) = -Vt ■ Q xy , (A20) 



Q' xz (z = 00) = - J - 2t ■ (Q xz - S voi sin Goo sin© c 
V z + K21 + K31 



(A21) 



e^ = -) = -V^ra' Q »" (A22) 

if,'(z = 00) = - J--ip. (A23) 

V K 

The Cauchy boundary condition at infinity has to be derived from variational calculus in 
order to be included into the relaxation. It is obtained from the bulk energy functional T 
and a boundary energy functional T^-, 

^ + U = ° > & = *>Q*)- ( A 24) 

To obtain this boundary energy functional we insert the Cauchy boundary conditions (A18)- 
(A23), previously derived from the linearized Euler-Lagrange equations, into the variational 
equation (A24). This immediately yields the six partial derivatives of the boundary energy 
functional, 





1 / 2 + k 21 + k 31 1 + k 2 i + hi 

2 \ 1 + k 21 + k 31 2 + k 21 + k 31 

1 / 2 + k 21 + k 31 1 + k 21 + k 31 
4 \ 1 + k 21 + k 31 2 + k 21 + k 31 

va + b + va-b va + b-va-b 
t/aTb -t/a^b t/aTb + ^a^b 

Qxx — SVoi (sin 2 — |) 

Qyy ~t~ 3 *~Vol 



(A25) 



(A26) 



U V 
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Qxx S vo \ (sin 2 |) 
^ Qyy 3 'S'voi 



(A27) 
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d&T y = - Q '*y = St-Q*y (A28) 



d 3~ i /i 

— - = -- (2 + k 21 + fc 31 ) = J- (2 + fc 21 + fc 31 ) t • (Q^ - S vo i sin Ooo sin 0^) (A29) 



_^ = _i ( 2 + ka + fai) Qi, = y ; < 2 + fei + *3i) * ' Q„* (A30) 

°-w = " kV/ (A31) 

= v 7 ^ 1 - ^ (A32) 

Now the boundary energy functional is obtained by integration, just in the same way a 
potential field is calculated from a given force field in classical mechanics. The final result, 
quoted in the main text (11), reads 

Foo = -U (Q xx 2 + Q yy 2 ) + - 3 (U + V) S vol (Q xx + Q yy ) (A33) 



-SVoi sin 2 Ooo • (U Q xx + V Q yy ) + V Q xx Q yy + -\/t-Q. 



2 

xy 



+ 2 V 2 ^ + + ^ t ' 2 ~ 2 Sy ° l Sil1 ' ^ + ^ 

2 
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FIG. 1: Profile of the smectic order parameter ip(z) for reduced temperature t = 0.0, at various 
intrinsic smectic wave numbers qo- Solid line: ip(z) at qo = 0.3; dashed line: ip(z) at qo = 0.6; 
dotted line: ip{ z ) at Qo = 0.8. 

FIG. 2: Profile of the nematic order parameters S(z) and T(z) for uniaxial and biaxial order, 
respectively, for reduced temperature t = 0.0, at various intrinsic smectic wave numbers qo. Solid 
line: S(z) at qo = 0.3; dashed line: S(z) at qo = 0.6; upper dotted line: S(z) at qo = 0.8; lower 
dotted line: T(z) at qo = 0.8. 

FIG. 4: Profile of the director tilt angle Q(z) for reduced temperature t = 1.0, at various intrinsic 
smectic wave numbers qo- Solid line: Q(z) at qo = 0.3; dashed line: 0(z) at qo = 0.6; dotted line: 
Q(z) at g = 0.8. 



FIG. 3: Profile of the director tilt angle Q(z) for reduced temperature t = 0.0, at various intrinsic 
smectic wave numbers qo. Solid line: ®(z) at go = 0.3; dashed line: G(z) at qo = 0.6; dotted line: 
0{z) at q = 0.8. 
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TABLE I: Tilt angle at infinity Gqo [degrees] in dependence on the reduced temperature t ( 
and the intrinsic smectic wave number go (columns). Both t and qo are in reduced units. 
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J. Stelzer and R. Bernhard : Figure 3 



^3 




distance from surface [reduced units] 



J. Stelzer and R. Bernhard : Figure 4 
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